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Abstract In the present paper we obtain an anisotropic analogue of Durgapal- 
Fuloria [1] perfect fluid solution. The methodology consists of contraction of anisotropic 
factor A by the help of both metric potentials e'“' and e^. Here we consider same 
as Durgapal-Fuloria [1] whereas e‘' is that given by Lake [2]. The field equations 
are solved by the change of dependent variable method. The solutions set math¬ 
ematically thus obtained are compared with the physical properties of some of 
the compact stars, strange star as well as white dwarf. It is observed that all the 
expected physical features are available related to stellar fluid distribution which 
clearly indicate validity of the model. 

Keywords General relativity, anisotropic fluid, compact stars 


1 Introduction 

Few decades ago a new analytic relativistic model was obtained by Durgapal and 
Fuloria [1] for superdense stars in the framework of Einstein’s General Theory of 
Relativity. They showed that the model in connection to neutron star stands all the 
tests of physical reality with the maximum mass 4.17 Mq and the surface redshift 
0.63. Very recently Gupta and Maurya [3] presented a class of charged analogues of 
superdense star model due to Durgapal and Fuloria [1] under the Einstein-Maxwell 
spacetimes. The members of this class have been shown to satisfy various physical 
conditions and exhibit features (i) with the maximum mass 3.2860 Mq and the 
radius 18.3990 km for a particular interval of the parameter 1 < K < 1.7300, and 
(ii) with the maximum mass 1.9672 Mq and the radius 15.9755 km for another 
interval of the parameter 1 < K < 1.1021. Later on a family of well behaved 
charged analogues of Durgapal and Fuloria [1] perfect fluid exact solution was also 
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obtained by Murad and Fatema [3] where they have studied the Crab pulsar with 
radius 13.21 km. 

In a similar way we have considered a generalization of Durgapal and Fuloria [1] 
with anisotropic fluid sphere such that pr ^ pt, where pr and pt respectively are 
radial and tangential pressures of fluid distribution. The present work is a sequel 
of the paper where we have developed a general algorithm in the form of 
metric potential u for all spherically symmetric charged anisotropic solutions in 
connection to compact stars. However, in the present study without considering 
any anisotropic function we can develop algorithm by the help of metric potentials 
only and here lies the beauty of the investigation. Another point we would like 
to add here that till now, as far as our knowledge is concerned, no alternative 
anisotropic analogue of Duragapal-Fuloria [1] solution is available in the literature. 

In connection to anisotropy we note that it was Ruderman who argued 
that the nuclear matter may have anisotropic features at least in certain very 
high density ranges (> 10^® gm/cm^) and thus the nuclear interaction can be 

treated under relativistic background. Later on Bowers and Liang specifically 
investigated the non-negligible effects of anisotropy on maximum equilibrium mass 
and surface redshift. In this regard several recently performed anisotropic compact 
star models may be consulted for further reference [8l[^ ll0l[mil21IT^I141ll51ll6| . We 
also note some special works with anisotropic aspect in the physical system like 
Globular Clusters, Galactic Bulges and Dark Halos in the Refs. mm- 

As a special feature of anisotropy we note that for small radial increase the 
anisotropic parameter increases. However, after reaching a maximum in the interior 
of the star it becomes a decreasing function of the radial distance as shown by Mak 
and Harko [I1[2Q]. Obviously at the centre of the fluid sphere the anisotropy is 
expected to vanish. 

We would like to mention that algorithm for perfect fluid and anisotropic un¬ 
charged fluid is already available in the literature msm- As for example, we 
note that in his work Lake plin] has considered an algorithm based on the choice 
of a single monotone function which generates all regular static spherically sym¬ 
metric perfect as well as anisotropic fluid solutions under the Einstein spacetimes. 
It is also observed that Herrera et al. [22] have extended the algorithm to the case 
of locally anisotropic fluids. Thus we opt for an algorithm to a more general case 
with anisotropic fluid distribution. However, in this context it is to note that in 
the Ref. |5] we developed an algorithm in the Einstein-Maxwell spacetimes. 

The outline of the present paper can be put as follows: in Sec. 2 the Einstein 
field equations for anisotropic stellar source are given whereas the general solutions 
are shown in Sec. 3 along with the necessary matching condition. In Sec. 4 we 
represent interesting features of the physical parameters which include density, 
pressure, stability, charge, anisotropy and redshift. As a special study we provide 
several data sheets in connection to compact stars. Sec. 5 is used as a platform for 
some discussions and conclusions. 


2 The Einstein field equations 

In this work we intend to study a static and spherically symmetric matter distri¬ 
bution whose interior metric is given in Schwarzschild coordinates, a;* = 






3 


EHlEa 


ds^ = — r^{d6^ + siri?6d(jP') + . 

The functions u and A satisfy the Einstein field equations, 

= R'j - \R9 3- 


( 1 ) 

( 2 ) 


where k = Stt is the Einstein constant with G = 1 = c in relativistic geometrized 
unit, G and c respectively being the Newtonian gravitational constant and velocity 
of photon in vacua. 

The matter within the star is assumed to be locally anisotropic fluid in nature 
and consequently T® j is the energy-momentum tensor of fluid distribution defined 

= [{P + PrWvj - ptS'-j + (pr - Pt)d''dj ], (3) 

where u® is the four-velocity as = 5*4, 0® is the unit space like vector 

in the direction of radial vector, 0® = is the energy density, pr is the 

pressure in direction of 0® (normal pressure) and pt is the pressure orthogonal to 
9i (transverse or tangential pressure). 

For the spherically symmetric metric o, the Einstein held equations may be 
expressed as the following system of ordinary differential equations [25] 




- kT‘^2 = -kT^3 = 


V 

~2 


\'v' ^ 

4 4 


Kpr, 


(4) 

v' 

2r 

e ^ = Kpt, 

(5) 

= i^P, 


(6) 


r 

where the prime denotes differential with respect to radial coordinate r. 
The pressure anisotropy condition for the system can be provided as 


A = n{pt -pr) = 


v/ / /2 t xf 

\ V V V — A 


-A -A , (1-e-^) 

e-e -h ^5 

r 


(7) 


Now let us consider the metric potentials [I] in the following forms: 

_A _ 7 - lOGr^ - gV^ 

® ~ 7+ 14Gr2 + 7G2r4 ’ 


( 8 ) 


v = (9) 

where G is a positive constant and ^ is a function which depends on radial coor¬ 
dinate r. The nature of plots for these quantities are shown in Fig. 1. 

The above Eq. (7) together with Eqs. (8) and (9) becomes 


'7- lOGr^ - gV®‘' 


'C^r^ + 19G^r®‘ - 21Cr‘^ - 7' 

tp 

'8C‘^r^{Cr^ +5) 

7(l-hGr2)2 

i) 

7r(l + Gr2)3 

7(l-hGr2)3 

(10) 
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r/R 


Fig. 1 Variation of the metric potentials with radial coordinate r/R are shown in the above 
figure. Here the legends are as follows: (i) e" is plotted with dotted line for Her X-1 and short- 
dashed line for white dwarf, (ii) is plotted with continuous line for Her X-1 and long-dashed 
line for white dwarf 



Fig. 2 Variation of the anisotropy factor with radial coordinate r/R are shown in figure. Here 
the legends are as follows: (i) A is plotted with short-dashed line for Her X-1, (ii) A is plotted 
with continuous line for white dwarf 


3 The solutions for the model 

Here our initial aim is to find out the pressure anisotropic function A, which is zero 
at the centre and monotonic increasing for suitable choices of ip. However, Lake 
[2] imposes condition that ip should be regular and monotonic increasing function 
of radial coordinate r. 

Let us therefore take the form of ip as follows: 

iP = {l-a + Cr^f, (11) 

where a > 0. 

Substituting the value of ip from Eq. (11) in Eq. (10), we get 

^ _ _8 Q[2-h (16 - g) Cr^ - 5Q - 2] 

7 (i+(7r2)3(l-a-^C'r2)2 


( 12 ) 
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For a > 0 and 0 < Cr^ < +8a+272-(i6--Rr) ^ pressure anisotropy is finite 

as well as positive everywhere as can be seen in Fig. 2. 

By inserting the above value of A in the Eq. (12), we get 


+ igC^r^ - 21C'r2 - 7 


+ 


1 


7 - 10Cr2 - C'2r4 




8C^r^(Cr^ + + (16 - a)Cr^ - (5a + 2)] 


t/)' 


(1 + Cr2) 


r(l + Cr^)(7 - 10C'r2 - C‘^r‘^) 

■>p =(ffl3) 


(l + C'r2)(l-a + C'r2)2 


Now our next task is to obtain the most general solution of the differential Eq. 
(13). Here we shall use the change of dependent variable method. We consider the 
differential equation of the form 


y" + p(r)y' + q(r)y = 0. 


(14) 


Let y = yi he the particular solution of the differential Eq. (14). Then y = yiU 
will be complete solution of the differential Eq. (14), where 


U = a\+h\ J exp — J (p(r) + 


dr, 


where ai and bi are arbitrary constants. 

Again let us consider here that = (1 —a+Cr^)^ = tpoir is a particular solution 
of Eq. (13). So, the most general solution of the differential Eq. (13) can be given 
by 


2\2 


Ip = (1—a+Cr^) 


B + A exp < — 


+ 


8Cr^ 


where A and B are arbitrary constants. 
After integrating it, we get 


C®r® + igC^r^ - 21Cr‘^ - 7 

r(l + C'r2)(7-10C'r2-C2r4) ' r(l - a + C'r2) 

(15) 


dr > dr 


i> = ipo 


where 


B-A- 


{V’qI + ^> 02(1 - a + Cr'^) + 1 pa 3 '(par}\/ (^q 5 “ 2(4 + a)(l - a + Cr'^) - pJar) 


(1 — a + Cr'^)^ 

= (1-a + Cr^f, 


+ W(r) 


(16) 

(17) 


W{r) = log 




V’aS - (4 + a)(l - a + Cr'^) + ^/ipab\/(P’a 3 - 2(4 + a)(l - a + Cr"^) - pjar) 

(1 - a + Cr^)p>a5 


and A and B are arbitrary constants with 
“ 3(16-8a-Q:^) ’ 


V’q2 


24-2a+a^ 
3(16-8a-a^)^ ’ 


V’aS 


288+80a-10a^ + a^ 
3(16-8a-a^)‘' ’ 


( 18 ) 
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i>aA 


1536-384a+48a^-2a^ 

3(16-8ct-a2)3 


V’aS = (16 — 8a — a^). 


Using Eqs. (8), (12) and (16) the expressions for energy-density and pressure 
read as 


Kp _ 8{9 + 2Cr^ + C‘^r'^) 
C7 “ 7(l-f Cr2)3 



r/R 


Fig. 3 Variation of the density with radial coordinate r/R are shown in figure. Here the 
legends are as follows: (i) p is plotted with continuous line for Her X-1 (ii) p is plotted with 
dashed line for white dwarf 




Fig. 4 Variation of the density with radial coordinate r/R are shown in figure. Here the 
legends are as follows: (i) pr is plotted with continuous line for Her X-1 and short-dashed line 
for white dwarf in left graph (left panel), (ii) pt is plotted with continuous line for Her X-1 
and short-dashed line for white dwarf in right graph (right panel) 


and 




















( 20 ) 


Kpr _ 4(7 - lOCr^ - C'V'^) rv>pr(l - a + Cr'^f + 2V> 

C ~ 7(1 + Cr2)2 [ V’(l-a + C'r2) 


8(C'r2 + 3) 
7(1 + Cr2)2’ 


where 


V’pr 


_ A{l + Cr^) _ 

(1 - a + Cr2)4y^(7 - 10C'r2 - Cr‘^)' 


( 21 ) 




Fig. 5 Variation of the density with radial coordinate r/R are shown in figure. Here the 
legends are as follows: (i) Pr/p is plotted with short-dashed line for Her X-1 and continuous 
line for white dwarf (left panel), (ii) pt/p is plotted with short-dashed line for Her X-1 and 
continuous line for white dwarf (right panel) 


In the Figs. 3-5 we have plotted the nature of the above physical quantities 
which show viable features of the present model. 


3.1 Matching condition 

The above system of equations is to be solved subject to the boundary condition 
that radial pressure pr = 0 at r = i? (where r = R is the outer boundary of 
the fluid sphere). It is clear that m(r = R) = M is a constant and, in fact, the 
interior metric (2.1) can be joined smoothly at the surface of spheres {r = R), to 
an exterior Schwarzschild metric whose mass is same as above i.e. m(r = R) = M 


The exterior spacetime of the star will be described by the Schwarzschild metric 
given by 

= dr'^-r^{de‘^+sin^ed(lP') + [l-‘^^dt^. ( 22 ) 


Continuity of the metric coefficients gtt, grr across the boundary surface r = 
R between the interior and the exterior regions of the star yields the following 
conditions: 


_ A(m 


2M 


-1 






















( 24 ) 




where 'i/){r = R) = ipj^. 

Equations (23) and (24) respectively give 

K,f - ^ \8CR^{3 + CR^y 
~ 2 [ 7(1 + C7i?2)2 J ’ 

^ _ V7 - 10C-R2 - C'^R'i 
~ V7{1 +CR^)i,^R{^-n{R)y 

The radial pressure pr is zero at the boundary (r = R) provides 

B_ (1 + CR'^)V7 - 10C'i?2 - C2i?2 

A ~ 2(1 - a + C'i?2)3[(i -a + CR‘^){3 + CR^) - (7 - IQCR^ - C^R^] ‘ 


where 


^afl = (l-a + C^^")^ 


n{R) = 


{^al + V’a2(l - a + CRy + i>a3'>PaR}\/i>a5 “ 2(4 + a)(l - O + (7772) - 

( 1-0 +( 7772)3 


+lT(i7), 


ip/'p'i _ '*/’a4 I '0a5 — (4; + a)(l — O + (777^) + y/lpa3\/ V’a5 “ 2(4 + Q;)(l — a + CR'^) — tpaR 

^ ® (1 - a + C772)V>a5 


4 Some physical features of the model 

4.1 Regularity at centre 


The density p and radial pressure pr and tangential pressure pt should be positive 
inside the star. The central density at centre for the present model is 

po = p(r = 0) = (31) 

The metric Eq. (22) implies that C = is positive hnite. 

Again, from Eq. (20), we obtain 

Pr{r = 0) ^ 4A _^ 

C V7(l-ayA=o 7’ ^ ^ 

where pr(r = 0) > 0. 

This immediately implies that 

B ^ V7 {'<pal + 1 pa 2 {l - o) + ^/'a3(l “ «)^}\/“ 2(4 + a)(l - a) - (1 - a)2 
A ^ 6(1-ay (1-0)3 

, -tpa 4 , I'V'aS - (4 + o)(l - o) + “ 2(4 + o)(l - o) - (1 - q)^] , 

VVa5 (-1 - a)'tpa5 
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4.2 Causality conditions 


Inside the fluid sphere the speed of sound should be less than the speed of light 


i.e. 0 < Vsr = \/^ < 1 and 0 < I4t = \/^ < 1- Therefore 


— / 


V sr 


(l+Cr^) 


4(C^r^ + lOCr^ 


7)(1 + Cr^)(^^i-V>^2 

8(C'2r4 + 2C'r2 + 25) 


- 8{Cr'^ + 5)' 
(34) 



r/R 


Fig. 6 Variation of the sound velocity with radial coordinate r/R are shown in figure. Here 
the legends are as follows: (i) Vr is plotted with dash line for Her X-1, (ii) Vr is plotted with 
marker continuous line for white dwarf (iii) Vt is plotted with continuous line for Her X-1 (iv) 
Vt is plotted with dotted line for white dwarf 



r/R 


Fig. 7 Variation of the sound velocity with radial coordinate r/R are shown in figure. Here 
the legends are as follows: (i) is plotted with dotted line for Her X-l, (ii) is plotted 
with dashed line for white dwarf (iii) is plotted with marker continuous line for Her X-1 
(iv) is plotted with continuous line for white dwarf 
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V\t = (l+Cr^) 

where 


4(c2r^ + lOCr^ - 7)(1 + Cr^){i^ri - V’r2 - i’rs) - 8{Cr^ + 5) - V>r4 


'^pr — 


8(C'2r4 + 2C'r2 + 25) 

^(l + Cr^) 


(1 - a + C'r2)4^(7- 10C'r2 - Cr^) 


Iprl = [ 


4{3 - Cr^){l - a + Cr^y^Ppr , 
(l-a + C'r2)2 ' + Cr2)(7-10C'r2 - C2r4)J 


+ 


'>Pr2 = [■; 


2K 


'IprS = 


(1 — a + Cr2) 
8(3-Cr^) p 


+ 


(1 — Q + Cr ) Ippr -12 




+ 


(l-a + C'r2)2^p 


(7-10Cr2-C2r4)(l + Cr2)^(l-Q + C'r2) ' i, 

(1 + Cr2)(l - a + Cr2)V>r5 - Cr'^i^rei^ - 3a + SCr^) 


^r4 = [4a' 


(1 + Cr2)(l — a + Cr^y 
= [12C'^r‘‘ + 4(16 - a)C'r^ - (10a + 4)], 


i/,^6 = [4(7^/ + 2(16 - a)C'r^ - (10a + 4)]. 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 


The physical quantities related to the above equations are plotted in Figs. 6 
and 7. 


4.3 Well behaved condition 

The velocity of sound is monotonically decreasing away from the centre and it 
is increasing with the increase of density i.e. ^(■^) < 0 or (qj^) > 0 and 

(F(^) < 0 or (^ 3 ^) > 0 for 0 < r < i?. In this context it is worth mentioning 
that the equation of state at ultra-high distribution has the property that the 
sound speed is decreasing outwards m as can be observed from Fig. 6. 


4.4 Energy conditions 

The anisotropic fluid sphere composed of strange matter will satisfy the null en¬ 
ergy condition (NEC), weak energy condition (WEC) and strong energy condition 
(SEC), if the following inequalities hold simultaneously at all points in the star: 
NEC: p> 0, 

WEC: p -|- Pr > 0, 

WEC: p-hpt > 0, 

SEC: p + pr + ‘2pt > 0. 

We have shown the energy conditions in Eig. 8 for Her X-1 under (i) and for 
white dwarf under (ii). 
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r/R 


Fig. 8 Variation of the energy conditions with radial coordinate r/R are shown in figure. 
Here the legends are as follows: (i) NEC is plotted with marker long-dashed line, WECr with 
marker short-dashed line and SEC with continuous line (for Her X-1), (ii) NEC is plotted 
with marker continuous line, WECr with long-dashed line and SEC with short-dashed line 
(for white dwarf) 


4.5 Stability conditions 
4-5.1 Case-1: 

In order to have an equilibrium configuration the matter must be stable against 
the collapse of local regions. This requires Le Chateliers principle, also known 
as local or microscopic stability condition, that the radial pressure pr must be a 
monotonically non-decreasing function of r such that > 0 Heintzmann 
and Hillebrandt also proposed that neutron star with anisotropic equation of 
state are stable for 7 > 4/3 as is observed from Fig. 9 and also shown in Tables 1 
and 2 of our model related to compact stars. 


4 . 5.2 Case-2: 


For physically acceptable model, one expects that the velocity of sound should be 
within the range 0 = = {dpi/dp) < 1 |3f)ll31| . We plot the radial and transverse 

velocity of sound in Fig. 7 and conclude that all parameters satisfy the inequalities 
0 = Vgr = (dpi/dp) < 1 and 0 = = (dpi/dp) < 1 everywhere inside the star 

models. Also 0 = Vg^ < 1 and 0 = Vgr < 1, therefore — 14^1 < 1- Now, to 
examine the stability of local anisotropic fluid distribution, we follow the cracking 
(also known as overturning) concept of Herrera |3U| which states that the region 
for which radial speed of sound is greater than the transverse speed of sound is a 
potentially stable region. 

For this we calculate the difference of velocities as follows: 


V^t - Vl = a 


(1 -h C'r2)(l -a + Cr'^)iPr5 - C'rVr6(5 - 3a + hCr^) 
2(1 -b C'r2)3(l -a + Cr2)3(c2j.4 2Cr2 -b 25) 


(43) 


where 


^^5 = [12C'^r'‘ + 4(16 - a)Cr'^ - (lOa + 4)], 


(44) 
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r/R 


Fig. 9 Variation of the adiabatic index with radial coordinate vjK are shown in figure. Here 
the legends are as follows: (i) 7 is plotted with continuous line for Her X-1 (ii) 7 is plotted 
with dash line for white dwarf 


= [AC'^r^ + 2(16 - a)Cr‘^ - (10a + 4)]. (45) 



r/R 


Fig. 10 Variation of the absolute value of square of sound velocity with radial coordinate r/R 
are shown in figure. Here the legends are as follows: (i) is plotted with short-dashed 

line for Her X-1, (ii) — Vgj.\ is plotted with continuous line for white dwarf 


It can be seen that — Vsr\ at the centre lies between 0 and 1 (see Fig. 
10). This implies that we must have 0 < < 1. Then a should satisfy the 

following condition: 0 < a < . 
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4.6 Generalized TOV equation 


The generalized Tolman-Oppenheimer-Volkoff (TOV) equation 


dr r 


( 46 ) 


where Mq = Mq (r) is the effective gravitational mass which can be given by 


i 9 / 

Mci'i') = 2 ^ ® ^ ^ ■ 


(47) 


Substituting the value of MQ{r) in Eq. (46), we get 

-\^'iP + Pr)+‘^{pt-Pr) = G. (48) 

Equation (48) basically describes the equilibrium condition for an anisotropic 
fluid subject to gravitational (Eg), hydrostatic (E/j) and anisotropic stress (Ea) 
which can, in a compact form, be expressed as 

Fg + Fh+Fa = 0, (49) 

where 

Fg = -\^'{p + Pr), (50) 


Fh = 


dpr 
dr ’ 


(51) 


2 

Fa = - {pt - Pr) . 
r 

The above forces can be expressed in the following explicit forms: 

^ 1 ,, , C^r \8{6-2Cr‘^)[i;pr{l-a + Cr^f+ 2i;] 

F, = - jP if + Vr) = ^ [ 7(i^cr»)3 - v(l-„ + Cp3) - 

4(7 - - lOGr^) f-iPpril - a + Cr'^f + 2V>y 

7(1 + Gr2)2 2 /’(l-a + C'r2) ) ’ 


(52) 


(53) 


Fh 


dpr 

dr 


C^r 

47r 


4(gV^ + lOCr^ 

7(1 + Gr2)2 


7) 


(' 0 rl 


^r2 - V’r-s) 


8 {Cr‘^ + 5) ■ 

7(1 + Gr2)3j ’ 

(54) 


Fa 


- (Pt - Pr) 
r 


C^r 


2a[(5a + 2) - (16 - a)Cr‘^ - 2C^r^] 
7(l + Gr2)3(l-a + Gr2)2 


(55) 


Variation of different forces and attainment of equilibrium has been drawn in 
Fig. 11. 
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Fig. 11 Variation of the forces with radial coordinate rjR are shown in figure. Here the 
legends are as follows: Fg is plotted with long-dashed line, with continuous line and Fa 
with short-dashed line 


4.7 Effective mass-radius relation and surface redshift 


Let us now turn our attention towards the effective mass to radius relationship. 
For static spherically symmetric perfect fluid star, Buchdahl [32] has proposed 
an absolute constraint on the maximally allowable mass-to-radius ratio (M/R) for 
isotropic fluid spheres as 2M/R < 8/9 (in the unit c = G = f). This basically states 
that for a given radius a static isotropic fluid sphere cannot be arbitrarily massive. 
However, for more generalized expression for mass-to-radius ratio one may look at 
the paper by Mak and Harko [9] . 

For the present compact star model, the effective mass is written as 


nR 

Mf.ff = 4 :Tv I pr^dr 

Jo 




SCR^iS + CR^) 
7(1-h 2Gi?2 G2i?4) 


(56) 


The compactness of the star is therefore can be given by 

_ M^ff _ 1 r 8CR‘^{3 + CR^) 

“ “ R ~ 2 7(1 + 2Gi?2 G2i?4) 


(57) 


Therefore, the surface redshift (Z) corresponding to the above compactness 
factor (u) is obtained as 


Z =[l- 2u]"^/^ - 1 = [1 


8CR^{3 + CR^) 
7{1 + 2CR^ + C‘^R^) 


(58) 


We have shown the variation of physical quantities related to Buchdahl’s mass- 
to-radius ratio (2M/R) for isotropic fluid spheres and also surface redshift are 
plotted in Figs. 12 and 13. 


5 Model parameters and comparison with some of the compact stars 

In this Section we prepare several data sheets for the model parameters in the fol¬ 
lowing Tables 1-3 and compare those with some of the compact stars, e.g. Strange 
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Fig. 12 Variation of the mass with radial coordinate r/R are shown in figure. Here the legends 
are as follows: (i) 2MJR is plotted with short-dashed line for Her X-1 and long-dashed line for 
its upper bound (left panel), (ii) 2M/R is plotted with short-dashed line for white dwarf and 
long-dashed line for its upper bound (right panel) 



r/R 


Fig. 13 Variation of the redshift index radial coordinate r/R are shown in figure. Here the 
legends are as follows: Z is plotted with continuous line for Her X-1 and long-dashed line for 
white dwarf 


star Her X-1 and White dwarf in Table 4. In our present investigation we propose a 
stable model with the parameters R = 12.5202 Km and M = 0.8882 Mq (for white 
dwarf) whereas R = 7.7214 Km and M = 0.8804 Mq (Her X-1) type compact star. 
The values of these data points have already been used for plotting graphs in the 
previous Sections 3 and 4 (See Figs. 1-13) in some way or others. 

Practically what we have done in the tables are as follows: In Tables 1-3 values 
of different physical parameters of Strange star Her X-1 and White dwarf have 
been provided. Under this data set then we calculate some physical parameters 
of compact star, say central density, surface density, central pressure etc in Table 
4. It can be observed that these data are quite satisfactory for the compact stars 
whether it is strange star with central density 1.0913 x 10^® gm/cmT^ or white 
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Table 1 Values of different physical parameters of White dwarf star for 
a = 0.10, CR^ = 0.068, M = 0.8882 M©, R = 12.5202 Km 


r 

Pr 

Pt 

P 

Vr 

Vt 

zd 

Z 

7 

0.0 

0.6386 

0.6386 

10.2857 

0.6135 

0.6034 

0.00000 

0.2036 

10.4949 

0.1 

0.6314 

0.6316 

10.2663 

0.6133 

0.6033 

0.00024 

0.2028 

10.5863 

0.2 

0.6096 

0.6106 

10.2084 

0.6127 

0.6031 

0.00093 

0.2003 

10.8718 

0.3 

0.5738 

0.5758 

10.1129 

0.6116 

0.6026 

0.0020 

0.1961 

11.3897 

0.4 

0.5245 

0.5278 

9.9814 

0.6100 

0.6019 

0.0034 

0.1903 

12.2191 

0.5 

0.4624 

0.4673 

9.8157 

0.6079 

0.6010 

0.0049 

0.1830 

13.5139 

0.6 

0.3885 

0.3949 

9.6185 

0.6053 

0.5997 

0.0064 

0.1741 

15.5930 

0.7 

0.3040 

0.3118 

9.3926 

0.6022 

0.5980 

0.0078 

0.1637 

19.2098 

0.8 

0.2101 

0.2190 

9.1412 

0.5985 

0.5959 

0.0089 

0.1519 

26.6355 

0.9 

0.1083 

0.1178 

8.8676 

0.5942 

0.5932 

0.0095 

0.1388 

49.2412 

1.0 

0.0000 

0.0096 

8.5754 

0.5893 

0.5899 

0.0096 

0.1244 

oo 


Table 2 Values of different physical parameters of Strange star Her X-1 for 
a = 0.11, = 0.1178, M = 0.8804 Mq, R = 7.7214 Km 

r 

Pr 

Pt 

P 

Vr 

Vt 

Z\ 

Z 

7 

0.0 

1.2135 

1.2135 

10.2857 

0.6730 

0.6624 

0.0000 

0.4010 

4.2917 

0.1 

1.1984 

1.1988 

10.2521 

0.6726 

0.6622 

0.0004 

0.3991 

4.3223 

0.2 

1.1533 

1.1551 

10.1523 

0.6713 

0.6617 

0.0018 

0.3933 

4.4180 

0.3 

1.0795 

1.0833 

9.9890 

0.6692 

0.6607 

0.0038 

0.3838 

4.5919 

0.4 

0.9790 

0.9851 

9.7665 

0.6662 

0.6592 

0.0061 

0.3707 

4.8710 

0.5 

0.8546 

0.8630 

9.4906 

0.6622 

0.6571 

0.0084 

0.3541 

5.3078 

0.6 

0.7096 

0.7197 

9.1681 

0.6572 

0.6541 

0.0101 

0.3342 

6.0116 

0.7 

0.5475 

0.5586 

8.8067 

0.6510 

0.6501 

0.0111 

0.3113 

7.2403 

0.8 

0.3725 

0.3835 

8.4143 

0.6436 

0.6450 

0.0110 

0.2856 

9.7716 

0.9 

0.1886 

0.1982 

7.9990 

0.6350 

0.6385 

0.0095 

0.2574 

17.4995 

1.0 

0.0000 

0.0068 

7.5686 

0.6249 

0.6305 

0.0068 

0.2271 

- 


Table 3 Values of the model parameters A, B, C and a for different compact stars 


Compact star 
candidates 

M 

(Me) 

R 

(Km) 

A 

B 

C 

a. 

White dwarf 

0.8882 

12.5202 

-2.1463 

0.5533 

4.3380 X 10-13 

0.10 

Her X-1 

0.8804 

7.7214 

-1.6255 

0.5301 

1.9758 X 10-13 

0.11 


Table 4 Energy densities, central pressure and Buchdahl condition 

for different compact star candidates for the above parameter values of Tables 1-3 


Compact star 

Central Density 

Surface density 

Central pressure 

Buchdahl condition 

candidates 

{gra/cra~‘^) 

(gm/cm~^) 

[dynej cm~‘^) 

{2M/R < 8/9) 

White dwarf 

2.3961 X lOii 

2.0 X lOii 

1.3392 X 103i 

0.1418 

Her X - 1 

1.0913 X 10i3 

0.8031 X 10i3 

1.1591 X 1033 

0.2280 


dwarf with central density 2.3961 x 10^^ gmjcm Likewise this feature of compact 
stars can be explored for some other physical parameters also. 
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6 Discussion and Conclusion 

In the present work we have investigated about an anisotropic analogue of Durgapal- 
Fuloria [I] and possibilities of interesting physical properties of the proposed model. 
As a necessary step we have contracted the anisotropic factor A by the help of 
both metric potentials e'^ and e^. However, is considered here same as Durgapal- 
Fuloria [1] whereas e'^ is that given by Lake [2]. 

The held equations are solved by the change of dependent variable method 
and under suitable boundary condition the interior metric ( 2 . 1 ) has been joined 
smoothly at the surface of spheres (r = R), to an exterior Schwarzschild metric 
whose mass is same as m(r = R) = M [26]. The solutions set thus obtained are 
correlated with the physical properties of some of the compact stars which in¬ 
clude strange star as well as white dwarf. It is observed that the model is viable 
in connection to several physical features which are quite interesting and accept¬ 
able as proposed by other researchers within the framework of General Theory of 
Relativity. 

As a detailed discussion we would like to put forward here that several ver¬ 
ification scheme of the model have been performed and extract expected results 
some of which are as follows: 

(1) Regularity at centre: The density p and radial pressure pr and tangential 
pressure pt should be positive inside the star. It is shown that the central density 
at centre is po = p{r = 0) = and pr{r = 0) > 0. This means that the density p 
as well as radial pressure pr and tangential pressure pt all are positive inside the 
star. 

(2) Causality conditions: It is shown that inside the fluid sphere the speed of 

sound is less than the speed of light i.e. 0 < I 4 r = < li 0 < 14 * = \J^ < 

1 . 

(3) Well behaved condition: The velocity of sound is monotonically decreasing 
away from the centre and it is increasing with the increase of density as can be 
observed from Fig. 6 . 

(4) Energy conditions: From Fig. 9 we observe that the anisotropic fluid sphere 
composed of strange matter satisfy the null energy condition (NEC), weak energy 
condition (WEC) and strong energy condition (SEC) simultaneously at all points 
in the star. 

(5) Stability conditions: Following Heintzmann and Hillebrandt [2^ we note 
that neutron star with anisotropic equation of state are stable for 7 > 4/3 as is 
observed in Tables 1 and 2 of our model. Also, it is expected that the velocity of 
sound should be within the range 0 = 17^ = [dpi/dp) < 1 |301I31| . The plots for 
the radial and transverse velocity of sound in Fig. 7 everywhere inside the star 
models. 

( 6 ) Generalized TOV equation: The generalized Tolman-Oppenheimer-Volkoff 
equation describes the equilibrium condition for the anisotropic fluid subject to 
gravitational (Fg), hydrostatic (F^) and anisotropic stress (Fa). Fig. 8 shows that 
the gravitational force is balanced by the joint action of hydrostatic and anisotropic 
forces to attain the required stability of the model. However, effect of anisotropic 
force is very less than the hydrostatic force. 

(7) Effective mass-radius relation and surface redshift: Eor static spherically 
symmetric perfect fluid star, the Buchdahl |32j absolute constraint on the maxi- 
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mally allowable mass-to-radius ratio (M/R) for isotropic fluid spheres as 2M/R < 
8/9 = 0.8888 is seen to be maintained in the present model as can be observed 
from the Table 4. 

In Sec. 5 we have made a comparative study by using model parameters and 
data of two of the compact stars which are, in general, very satisfactory as com¬ 
pared to the observational results. However, at this point we would like to comment 
that the sample data used for verifying the present model are to be increased to 
obtain more satisfactory and exhaustive features in the realm of physical reality. 
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